Dynamics of rigid and flexible extended bodies in viscous films and membranes 
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We study the dynamics of extended rod-like bodies in (or associated with) membranes and films. 
We demonstrate a striking difference between the mobilities in films and bulk fluids, even when the 
dissipation is dominated by the fluid stress: for large inclusions we find that rotation and motion 
perpendicular to the rod axis exhibits purely local drag, in which the drag coefficient is algebraic in 
the rod dimensions. We also study the dynamics of the internal modes of a semiflexible inclusion 
and find two dynamical regimes in the relaxation spectrum. 
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The mobility of inclusions in membranes is a funda- 
mental physical parameter controlling a number of cel- 
lular processes. Since inclusions such as proteins Q, Q 
or "rafts" 0, 0] can in many cases be large compared 
with the constituent lipids of the membrane, they can 
be viewed as macroscopic objects moving in a contin- 
uum fluid environment. The dynamics of these objects 
in thin films has already been shown to be surprisingly 
subtle, leading to confusion in the early literature on 
protein diffusion in cell membranes This contro- 

versy was clarified by Saffman |(| , who noted that motion 
of objects in or attached to a (two dimensional) mem- 
brane is never strictly two-dimensional since in-plane mo- 
mentum induces flows in the surrounding bulk (three- 
dimensional) fluids due to a viscous coupling of the in- 
terface/membrane to the surrounding fluids. That cou- 
pling introduces a new length-scale £q = r) m /r)t deter- 
mined by the ratio of membrane and fluid viscosities, 
?7m,f- This length determines the degree to which the 
dissipation is predominantly two- or three-dimensional 

n ie ii ii s m As a result, the drag coefficient on 

a small object in a membrane is a nonlinear function 
of both its size and the membrane viscosity. For exam- 
ple, the diffusion coefficients of small proteins (linear size 
< £o) in membranes have a weak (logarithmic) depen- 
dence on their size. 

Here, we examine the motion of rod-like inclusions in 
viscous films and membranes, as a representative exam- 
ple of a general scheme for the calculation of mobilities 
of arbitrary extended bodies, which we also present. The 
generalization of this problem to motion in viscoelastic 
films is straight-forward. Our work is motivated in part 
by prospect of computing mobilities of proteins and lipid 
rafts in lipid bilayers Q, [3] , as well as by recent exper- 
iments that have demonstrated the possibility of mak- 
ing quantitative rheological measurements of viscoelas- 
tic films resembling cellular structures such as the actin 
cortex |tj using tracer particle fluctuations (membrane 
microrheology) . Given its flexibility, this formalism can 
be adapted to the computation of the mobility of highly 
complex and irregular interfacial objects such as fractal 
aggregates [lOj. Finally, the driven motion of rods in vis- 
cous/ viscoelastic films has also been used to determine 



rheological properties, e.g., of monolayers |ll|,for which 
the approach we develop here is important. We also an- 
alyze the undulatory motion of semiflexible polymers at 
a membrane surface. 

Our calculation of the hydrodynamic drag on a rigid 
rod of length L, gives two principal results: (i) for small 
objects (specifically, L <C £o), the drag coefficients be- 
come independent of both the rod orientation and aspect 
ratio; and (ii) for larger rods of high aspect ratio, (± be- 
comes purely linear in the rod length L — i.e., the drag 
becomes purely local. In contrast, we find that the well- 
established three-dimensional result C\\ = 2Trr] / ln(AL / a) 
applies for parallel motion in the film, provided that 
L ^> £q. Here, however, the effective rod radius be- 
comes £q rather than the physical radius a, when a <C £q- 
Closely related to (ii), we find that the rotational drag 
(cquivalently diffusion constant) depends purely alge- 
braically on the rod length. 

We also calculate the relaxation spectrum of undula- 
tion modes of wavelength q^ 1 of a semiflexible inclusion 
of rigidity k and find a dynamic relaxation time r(q) 
which crosses over from r ~ q~ 3 on modes with wave- 
length less than £q to r ~ q~ 4 for those longer than £q. 
The latter result represents purely local drag, with no 
additional logarithmic dependence in r(q), as there is for 
bending fluctuations in a bulk fluid. 

The starting point of our analysis will be the re- 
sponse of the two dimensional membrane fluid to an 
in-plane force distribution. This is calculated by solv- 
ing the coupled equations for in-plane and out of plane 
fluid motions taking into account incompressibility |l3| 
of both the bulk and the membrane leading to an ex- 
pression for the in-plane velocity v a (x) at position x re- 
sulting from a point force fp (x') at another point x': 
v a (x) = a a p (x — x') fp (x') . The response function 
a a p (x) is given in closed form as 

a af3 (x) = a\\ (\x\)x a xp + a± (|x|) [5 a p - x a xp] , (1) 

where the scalar functions a\\ , a± of the distance between 
the point of the force application and the measurement 
of the velocity field are given by: 

Hi(z) 2 Y (z) + Y 2 {z) 



-4r] m iu!a\i(x, uj) 
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-4r] m iu)aj_(x,u)) = 

zH (z)-Hi(z) 2 Y Q (z)-Y 2 (z) 

z + ttz* 2 {6) 

where the are Struve functions 0], and the Y v are 
Bessel functions of the second kind. Here, z = |x| /£q is 
the distance between the point of force application and 
the membrane velocity response measured in the flat film 
in units of £q. 

To parameterize this geometry of the rod of length L 
and circular cross-section a, we define the dimensionless 
aspect ratio p = L/a. There are three independent drag 
coefficients to determine in the problem. The in-plane 
translational mobility tensor, p a p (and consequently, its 
inverse, the drag ( a = p,~^) is defined by 

„,rod ,, prod t A \ 

V a = VapFp , (4) 

where v T a od is the a th component of the velocity of the 
rod and F^ od is the /3 th component of the total force 
applied to the rod (at its center) and a, f3 € {1,2}. In- 
plane rotational symmetry and inversion of the rod imply 
that the mobility tensor has the form: /x^ = /iiinjfij + 
A*_l {Sij — fiifij) . Here fj,j_ and [in are the mobility of the 
rod dragged perpendicular to and parallel to its long axis 
(oriented along the n direction) respectively. In addition 
to these two independent translational mobilities, there 
is also one rotational mobility, /x ro t linking the angular 
velocity of the rod to the torque applied to that rod about 
its center of inversion symmetry. 

We solve the problem by two complementary methods 
useful for the regimes of small and large aspect ratios 
respectively. For small p, we use a 2-d analogue of the 
Kirkwood approximation |15| to model the continuous 
rod by a series of discs subject to point forces at their 
centers. This method becomes rather cumbersome when 
p ;§> 1, but here we may proceed by a second approxima- 
tion that assumes the rod to be infinitely thin. In both 
cases, we restrict our attention to the limit a <C £q. Be- 
low, we illustrate both methods by a calculation of the 
transverse drag coefficient with the understanding that 
the longitudinal and rotational drag proceed analogously. 

The linearity of the underlying low-Reynolds number 
hydrodynamics allows to use superposition. Specifically, 
we replace the rod of length L and cross-sectional radius 
a by a set of N + 1 disks of radius a and interdisk sep- 
aration b chosen so that the total length of the rod is 
preserved, i.e. L — Nb + 2a. We choose the number of 
disks to be maximal consistent with a given aspect ratio 
and the non-interpenetrability of the disks. 

Our strategy for computing the drag on the rod in- 
volves setting the rod in uniform motion with unit veloc- 
ity by imposing some set of forces fW, i = 1, . . . ,N + 1 
on the N + l disks making up the rod. We can use the re- 
sponse function to compute the velocity field for a given 



collection of point forces vfi ± — a \\±fn ±- How- 

ever, we must demand that all the disks have the same ve- 
locity and thereby determine the forces applied to them. 
To enforce this constraint, we invert the matrix ■ 

The drag coefficient is then, = Yn~jli ("I^l) 

For large N, i.e. high aspect ratio rods, the matrix 
inversion becomes difficult. To study that limit, one can 
assume an infinitely thin rod. Then, the velocity field 
at the point x due to a continuous distribution of force 
densities along the rod, f(xx) that lie along the x-axis 
from x = — L/2 to x — L/2 takes the form 

,L/2 

v a (x) = / a a/3 (x - px) fp (px) dp . (5) 

J -L/2 

As before, we impose a unit velocity field on the rod and 
determine the force density required to effect this result, 
f (x). The inversion of Eq.[S] proceeds by first expanding 
the linear force density in Legendre polynomials P n (x) 
writing 

N 

/(a;) = ^c 2 „P 2 „(2x/L) (6) 

71 = 

where the coefficients are as yet unknown. Since the Leg- 
endre polynomials form a complete set on the interval 
— I to +1, any physical force density can be expressed 
as in Eq. provided N be taken to infinity. In practice, 
we find excellent numerical results even when truncating 
this sum to just the first five terms, while taking into ac- 
count the symmetry of the force distribution about the 
center of the rod. 

We now impose the unit velocity condition at a fi- 
nite set of points < pi < L/2 along the rod since 
if we truncate the Legendre function expansion of the 
force density at N, we can impose the velocity condition 
at a maximum of N points without creating an over- 
determined system of equations. Thus we require that 
v (pix) = 1, fori = 1, . . . ,N. Now finding the force dis- 
tribution along the rod requires only the inversion of an 
N x N matrix, Afij whose components are defined by 

,L/2 

Afij = I a (xp t - zx) Pj(2z/L) dz. (7) 

J-L/2 

Finally, the total force on the rod is found by recon- 
structing the force density from its Legendre polynomial 
expansion and integrating the resulting expression over 
the entire rod. That force density is given by Eq.[S] where 
the coefficients, Ck are determined using the rod's unit ve- 
locity condition enforced at a discrete set of points along 
with the inverse of the matrix defined in Eq. [7] Thus we 
find Ck — X^i-^fci 1 - D ue to the orthonormality of the 
Legendre polynomials, the total force is given entirely 
by the coefficient of the zeroth Legendre polynomial, and 
the total torque in the case of rotations is given by the 
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coefficient of the first odd Legendre polynomial, c\. Nu- 
merically, we find a variation of less than a percent in 
the drag coefficient in all cases, so long as N > 2. Here, 
we report our results for N — 5. To check for consis- 
tency, we compared our results from this method and the 
Kirkwood approximation; for thin rods, we find excellent 
agreement. 
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FIG. 1: Parallel (||) and perpendicular _L drag coefficients for 
thin rods (i.e. high aspect ratio) of various lengths. For short 
rods both the parallel and perpendicular drag coefficients ap- 
proach the Safrman-Delbriick result. The divergence between 
these two drag coefficients at long lengths reflects the absence 
of the hydrodynamic cooperativity in perpendicular case. 

The analogous calculation can be made for the drag on 
the rod moving in a direction perpendicular to its long 
axis. The two drag coefficients (parallel and perpendic- 
ular) converge rapidly for rods of length L < £q, while 
their difference grows monotonically with increasing rod 
length for L > Iq, This results stands in contrast to the 
case of motion in bulk fluids, where the two drag coeffi- 
cients differ by a constant multiple. 

We plot the parallel and perpendicular drag coefficients 
for thin rods as a function of their length, L/£o in reduced 
units in Figure ^ The equality of the two drag coeffi- 
cients when all dimensions of the rod are small compared 
to £q reflects the effective shape and size independence of 
the Saffman-Delbriick result. For rods longer than this 
crossover length, these two drag coefficient diverge from 
each other while they both become more strongly length 
dependent. 

To understand this peculiarity of interfacial drag, we 
note two points. First, for long rods (L > Iq), the parallel 
drag in the film is essentially unchanged from the bulk, 
three-dimensional drag, in that C|| = \ n (Qi^L/i ) ' wnere 
the prefactor in the logarithm has been determined to 
within 1%. Comparing with the result for drag of a rod 
in a bulk fluid, we see that the effective radius of the rod 
is now of order £o (for a <C £q). Given that £q corresponds 
to a length scale over which interfacial momentum den- 
sities flow into the bulk fluid due to the viscous coupling 
between the two, the system does not effectively resolve 
length scales smaller than £q so the small dimension of 
the rod is replaced by this length. At large length scales 
(compared to £q) the fluid velocity field around the rod 
in parallel drag is the same as for rod motion in the bulk 
fluid, i.e. in both cases there should have been no flow 
from the plane of the interface into the subphase. 



The case of perpendicular motion of long rods, on the 
other hand, is qualitatively very different. Here we find 
C_l = 2ttt]L. For this case, the fluid flow field for a rod in 
bulk is inconsistent with the flow restrictions imposed by 
the presence of the interface. In 3d, there would be a non- 
vanishing two-dimensional divergence of velocity field re- 
stricted to the plane of motion of the rod. Now, the 
in-plane incompressibility of the interface requires that 
the fluid velocity field extend over distances comparable 
to the rod's largest dimension L. The standard hydrody- 
namic coupling of portions of the rod, which gives rise to 
the logarithm in the drag is not present resulting in a drag 
coefficient that is purely linear in rod length. In other 
words, the drag is effectively local or "free-draining" in 
character. In three dimensions, in contrast, there is a 
length-independent ratio of two between the parallel and 
perpendicular drag coefficients because neither is free- 
draining: C,? d = and c! d = 2$ d . 
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FIG. 2: The rotational drag coefficient of a rod of infinite 
aspect ratio plotted versus the length of the rod. 
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FIG. 3: The decay rate vs wavevector of transverse fluctu- 
ations of a semiflexible rod embedded in the interface. The 
decay rate scaled by q 4 is wavevector-independent for purely 
local drag (small q regime). Nonlocal hydrodynamic interac- 
tions al ong the rod reduce the decay rate in a manner analo- 
gous to |lal in the larger q regime. 

The calculation of the rotational drag coefficient pro- 
ceeds analogously to those of the perpendicular and par- 
allel drag coefficients. We plot the rotational drag coef- 
ficient divided by L 2 for a rod of infinite aspect ratio as 
a function of the reduced length in figure [3 The essen- 
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tial feature of this plot is that rotational drag coefficient 
scales as L 2 for rods smaller than £ and then as L 3 for 
rods longer than this natural length. Thus, we find purely 
algebraic behavior in both limits. 

Finally, to study the internal mode dynamics of a semi- 
flexible inclusion, we consider the small transverse fluc- 
tuations of a an infinitely thin, almost rod-like filament 
oriented along the x axis parameterized by (x, r(x, t)). 
With no slip boundary conditions on the inclusion we 
have a dynamical equation for r(x) to linear order, 

d t r(x, t)= dx'a(x - x', 0) • (0, -Kd*r) + r){x, t) (8) 

where the thermal fluctuations rj{x, 1) are chosen to sat- 
isfy the Fluctuation-Dissipation Theorem. Expanding in 
modes of wavenumber q, we obtain a relaxation spectrum 
of transverse correlations given by 



, , . . „ 2k B T , , , S1 
(r(q, t)f(q, 0)> = — cxp [-t/r(q)] 



nq 1 



(9) 



2nq 4 



dx 



where t 1 (q) - , 

^ ' Tim J q 2n X 2 {x 2 _ q 2y/2( x + e -i) 

We have obtained a complicated but closed form expres- 
sion for r(q) that is plotted in figure^ ft is simple in the 
two limits of q£ <C 1, qto 3> 1 where T~ 1 (q <C ^q" 1 ) — 

— — and r _1 (q 3> i^ 1 ) ~ — — respectively. At short 

scales we find the dimensionally reduced analogue of the 
hydrodynamic relaxation spectrum of fluid membranes ( 
D = 2 dimensional manifold embedded in D + 1 dimen- 
sions) 0| (we consider D = 1). Here though, the 'long- 
range' hydrodynamic coupling crosses over to a purely 
local friction on the longest length scales • 



In summary, we have calculated the translational and 
rotational hydrodynamic drag on a rod moving at low 
Reynolds number in a viscous film coupled to viscous 
sub- (and/or) superphase. When the dimensions of the 
rod are small (<SC £q), the dissipation is governed primar- 
ily by the film, and the drag is insensitive to orientation 
and aspect ratio, and only weakly (logarithmically) de- 
pendence on size 6]. We find, surprisingly, that in the 
limit of a long rod, the drag for motion perpendicular to 
the rod axis exhibits purely local drag per unit length. 
Thus, in contrast with motion in 3d fluids, there are 
no long-range hydrodynamic effects for long rods, even 
though the dissipation occurs entirely in the fluid. A 
similar observation applies to rotational motion, and for 
the relaxation spectrum of the long wavelength modes 
of a fluctuating filament or semi-flexible polymer in a 
viscous film. This is because, although the dissipation 
is governed primarily by the fluid, the film (along with 
its assumed incompressibility and no-slip conditions) im- 
poses a very different boundary condition on the flow 
from what we would have in a bulk fluid alone. This 
added condition not only increases the drag for perpen- 
dicular motion relative to that without the film present, 
but results in purely local, or free-draining drag dynam- 
ics. In contrast, since the new boundary conditions are 
consistent with the standard flow field for parallel mo- 
tion, we find quantitative agreement with the parallel 
mobility in a bulk fluid when the film's viscosity becomes 
irrelevant (small £q). Finally, we note that the methods 
developed here constitute a highly adaptable framework 
to compute the mobility of arbitrary extended, irregu- 
larly shaped objects embedded in a viscous membrane or 
interface. 
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